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ABSTRACT 


A  number  of  tests  and  estimates  for  mean  life  and 

other  parameters  derived  under  the  exponential  distribution 

assumption  are  studied  under  the  alternative  condition  that 

the  distribution  has  an  increasing  (decreasing)  failure  rate. 

The  tests  and  estimates  considered  are,  for  the  most  part, 

based  on  censored  and  truncated  samples.  It  is  shown  that 

the  usual  acceptance  sampling  procedures  based  on  mean  life 

generally  favor  the  producer  (consumer)  in  the  IFR  (DFR) 

case.  However,  acceptance  sampling  procedures  based  on  the 
t  h 

q  quantile  tend  to  favor  the  consumer  in  the  IFR  case  when 
q  is  small.  The  usual  estimates  for  the  mean  based  on  the 
exponential  assumption  are  positively  (negatively)  biased 
when  the  distribution  is  IFR  (DFR).  Many  of  these  results 
hold  under  the  weaker  assumption  that  the  distribution  has 
an  increasing  (decreasing)  failure  rate  on  the  average. 
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1.  INTRODUCTION 


In  a  fundamental  paper  in  the  literature  of  life  testing,  Epstein 
and  Sobel  (1953)  introduce  life  test  procedures  based  on  the  exponential 
distribution.  (See  also  Epstein  (1960a).)  These  procedures  have  been 
codified  in  a  Department  of  Defense  handbook  H108  (1960)  and  are  now 
widely  employed.  If  the  exponential  procedures  are  based  on  a  complete 
sample,  the  central  limit  law  assures  "robustness’1  for  large  sample 
sizes;  i.e.,  the  risk  function  is  not  sensitive  to  departures  from  the 
exponential  distribution  assumption.  However,  the  robustness  of 
exponential  procedures  based  on  censored  or  truncated  samples  has  been 
questioned.  Zelen  and  Dannemiller  (1961)  show  by  using  Weibull  distribution 
alternatives  that  these  procedures  may  not  be  robust  in  testing  for  mean 
life.  In  particular,  they  show  that  the  use  of  these  procedures  may  result 
in  substantially  increasing  the  probability  of  accepting  items  having  poor 
mean  lives.  Antelman  and  Savage  in  unpublished  work  point  out  that  these 
exponential  procedures  used  in  testing  hypotheses  concerning  the  median 
seem  to  be  robust.  This  is  not  the  case  in  testing  for  low  or  high 
percentiles.  Although  the  exponential  statistics  used  In  hypothesis  testing 
may  suffer  in  certain  applications  from  lack  of  robustness,  they  do  yield 
conservative  tolerance  and  confidence  limits  for  distributions  with 
increasing  failure  rate  in  certain  cases.  This  is  suggested  by  the  results 
of  Zelen  and  Dannemiller  and  developed  in  detail  in  Barlow  and 
Proschan  (1966b)  (hereafter  referred  to  as  BP  (1966b).) 
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Since  the  usual  tests  and  estimates  based  on  the  exponential 
assumption  possess  an  undeniable  computational  simplicity  and  elegance, 
perhaps  a  further  investigation  of  their  properties  when  the  true  dis¬ 
tribution  is  other  than  exponential  is  warranted.  This  paper 
interprets,  in  a  life  testing  context,  theoretical  results  concerning 
linear  combinations  of  order  statistics  contained  in  Barlow  and 
Proschan  (1966a)  (hereafter  referred  to  as  BP  (1966a)).  We  consider 
statistical  procedures  based  on  an  exponential  distribution  for  a 
variety  of  life  testing  plans.  In  each  case — censored  sampling, 
censored  sampling  with  replacement,  or  truncation — the  usual  estimate 
for  the  mean  life  is  biased  high  (low)  if  the  true  distribution  has  an 
increasing  (decreasing)  failure  rate.  Properties  of  distributions  with 
monotone  failure  rate  are  exploited  to  obtain  bounds  on  the  risk  function 
when  the  usual  decision  rules  based  on  an  exponential  assumption  are  used. 

Mathematical  Preliminaries.  Let  X  denote  a  random  variable  with 
right  continuous  distribution  F  such  that  F(0  )  =  0.  If  F  has  density 
f,  then  r(t)  *  f(t)/F(t)  is  known  as  the  failure  rate, where  F  =  1-F. 
Note  that  r(t)  =  -  ~  log[F(t)]  when  a  density  exists.  For  this  reason, 
we  say  that  F  is  IFR  (DFR)  for  increasing  (decreasing)  failure  rate  if 
log[F(t)]  is  concave  where  finite  [convex  on  [0,°°)).  Note  that  any 
IFR  (DFR)  distribution  with  specified  mean  can  be  expressed  as  the  limit 
of  continuous  IFR  (DFR)  distributions  with  the  same  mean.  Hence  for  many 
of  our  results  it  is  sufficient  to  confine  attention  to  continuous  IFR 
(DFR)  distributions  and  we  do  this  in  giving  proofs. 
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A  weaker  restriction  than  monotonicity  of  the  failure  rate  is 

also  considered  in  some  cases;  namely,  monotonicity  of  the  failure  rate 
1 

average,  —  /  r(u)du.  This  motivates  the  definition  of  IFRA  (DFRA) 
t 

distributions.  F  is  said  to  have  an  increasing  (decreasing)  failure 
rate  average}  i.e.,  to  be  IFRA  (DFRA)  if  -log  F(x)/x  is  increasing 
(decreasing).  Clearly,  if  F  is  IFR  (DFR) ,  then  it  is  IFRA  (DFRA). 

The  virtues  of  the  IFRA  class  of  distributions  are  discussed  in 
Birnbaum,  Esary  and  Marshall  (1966). 

Properties  of  distributions  with  monotone  failure  rate  are  discussed 
in  Barlow,  Marshall  and  Proschan  (1963)  and  in  Barlow  and  Proschan  (1965), 
Chapter  2.  Tables  of  bounds  for  these  distributions  are  given  in  Barlow 
and  Marshall  (1965).  Additional  bounds  are  given  in  Barlow  and  Marshall 
(1966).  Distribution  free  life  test  sampling  plans  based  on  these  bounds 
are  discussed  in  Barlow  and  Gupta  (1966). 

Unless  otherwise  indicated,  we  denote  ordered  observations  from  a 

random  sample  of  size  n  based  on  a  random  variable  X  by  X,  <•••<  X 

l,n  —  ~  n,n 

and  define  Xq  n  =  0.  The  second  subscript  is  dropped  when  there  is  no 
danger  of  confusion.  The  term  increasing  (decreasing)  is  used  for 
nondecreasing  (nonincreasing).  It  will  be  convenient  to  let  G  denote 
the  exponential  distribution  with  mean,  6,  i.e. 


1 


x  <  0 


G(x)  =  { 


x 

0 


e 


x  >  0. 


4 


We  shall  consider  distributions  F  that  are  IFR  (DFR)  or 
IFRA  (DFRA)  with  an  unknown  scale  parameter  of  interest,  such  as 
the  mean  or  a  specified  quantile.  We  shall  compare  operating 
characteristic  (OC)  curves  for  the  F  distribution  with  corresponding 
OC  curves  for  the  exponential  distribution  as  the  scale  parameter 


varies. 
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2.  CENSORED  SAMPLING  TO  ESTABLISH  A  SPECIFIED  MEAN  LIFE 

Assume  n  items  are  placed  on  test  simultaneously  and  testing  is 
discontinued  after  the  (r  <  n)  failure.  There  is  no  replacement 
of  failed  items.  If  the  distribution  is  exponential  with  density 


then 


g(x;6)  =  \ 


6  6 


x 

0 


x  >_  0 

x  <  0, 


0 

rn 


r 

£  xi  +  (n-r)xr 

r 


(1  <  r  <  n) , 


the  total  time  on  test  divided  by  the  number  of  failures  observed,  is 
the  maximum  likelihood  ai*.d  minimum  variance  unbiased  estimate  for  0 
based  on  the  first  r  order  statistics  (Epstein  and  Sobel  (1953)). 

The  following  decision  rule  based  on  §rn  is  a  uniformly  most  powerful 
test  of  the  hypothesis  0  =  0q  against  0  <  0^  with  Type  I  error  equal 
to  a. 

9  — >  accept 

x„(2r) 

c<eo)’eo“V"  <2a> 

— ► reject 


a.  £■'?*. 


.  MRMNTti  *'  :  I 
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2 

where  xa(2r)  denotes  the  100a  per  cent  point  of  a  chi-square 
distribution  with  2r  degrees  of  freedom  and  0q  is  the  mean  life 
goal.  Zelen  and  Dannemiller  (1961)  show  that  if  the  distribution  is 
actually  Weibull  with  increasing  failure  rate,  then  with  n  ■  28, 
r  *  14,  9q  *  1000  hours,  and  a  *  .10,  there  is  a  high  probability 
of  accepting  items  with  a  mean  life  considerably  less  than  0^  ■  1000 
hours.  This  phenomenon  is  true  in  general  when  the  distribution  is 
actually  IFR,  as  stated  more  precisely  in 

THEOREM  2.1:  If  F  is  IFR  (DFR)  with  mean  0  and  G  is 
exponential  with  mean  0,  then 


V®™  l  c(e0) I e}  i  (±)  pG{ern  ic(0o>le}  (2-2) 

for  0  >_  eQx^(2r)/2(n-r+l) . 

Inequality  (2.2)  follows  from  the  result 


r 

(n-i+1)  (Xi-Xi-1)  >x)  >  Pg(£  (n-i+1)  (Yi-Yi_1)  >_  x) 


for  x  (n-r+l)0  and  F  IFR  which  is  proved  in  BP  (1966a), 

Corollary  5.5  and  Theorem  4.6. 

For  small  r,  the  coefficient  of  0Q  is  less  than  1  (e.g.,  if 
a  £e  1  and  r  <_  (n+l)/2,  then  [x^(2r)/2(n-r+l) ]  <_  1);  see  BP  (1966b). 
For  the  Zelen-Dannemiller  example  we  can  assert  that  the  probability  of 
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*?10(28) 

accepting  items  with  true  mean  lives  as  low  as  - — -  100%  =  63% 

of  the  mean  life  goal  is  greater  when  the  true  distribution  is  IFR 
than  when  it  is  exponential.  Hence  we  see  that  the  producer's  risk 
is  controlled  against  IFR  alternatives,  but  not  the  consumer's  risk. 

If  we  keep  r  fixed  and  increase  n,  inequality  (2.2)  is 
strengthened . 

THEOREM  2.2:  If  F  is  IFR  (DFR)  with  mean  6,  then 


PF{%  VVn+liC(Vle)  (2l3) 

for  all  0  >  0. 

Inequality  (2.3)  is  proved  in  BP  (1966a),  Theorem  5.4.  It  follows 

A 

from  this  result  that  0  as  an  estimate  for  0  is  positively 
(negatively)  biased  when  F  is  IFR  (DFR). 

COROLLARY  2.3:  If  F  is  IFR  with  mean  0,  then 


(2.4) 


If  F  is  DFR  with  mean  0,  then 


0  <  E_,[0  ,.]  <  E_,0  <  0 

—  F  r,n+l  —  F  rn  — 


(2.5) 


for  1  <  r  <  n. 
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All  Inequalities  are  sharp. 

Proof .  The  upper  bound  In  (2.4)  follows  from  the  fact  that 
for  every  sample  realization 

r  n 

£x.  +  (n-r)X  l£x  . 

1  1 

Since  equality  Is  attained  with  the  distribution  degenerate  at  6 
(which  is  the  limit  of  IFR  distributions)  the  bound  is  sharp.  The 
other  inequalities  in  (2.4)  follow  from  Theorem  2.2. 

The  lower  bound  in  (2.5)  is  sharp  as  shown  in  Section  6  of 
BP  (1966a).  The  other  inequalities  follow  from  Theorem  2.2.  |  | 

Actually  it  can  be  shown  under  the  weaker  condition  of  F 
IFRA  (DFRA)  that  E0rn  is  increasing  (decreasing)  in  r  for  fixed 
n  and  hence  E6rn  ^  (f)  6.  (See  Theorem  5.3,  BP  (1966a).) 
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3.  CENSORED  SAMPLING  TO  ESTABLISH  A  SPECIFIED  QUANTILE  LIFE 

Suppose  now  that  we  wish  to  establish  that  the  quantile, 

£  ,  exceeds  some  quantile  life  goal,  say  £°,  with  a  Type  I  error 
<_  a.  We  consider  a  censored  sampling  plan,  as  in  Section  2.  If  the 
distribution  is  exponential,  the  decision  rule  is,  using  (2.1)  and 
the  relationship  between  and  0, 


Oi?(2r) 

q  a 

-  2r  5.n(l-q) 


0  1  C’(£0)  = 

rn  \  sq 


accept 


reject. 


(3.1) 


See  Epstein  (1960c).  In  reliability  situations  we  will  usually  be 
interested  in  low  percentiles;  i.e.  small  values  of  q.  For  such 
values  of  q,  the  consumer  rather  than  the  producer  will  be  protected 
under  this  decision  rule,  if  the  distribution  is  actually  IFR.  More 
precisely,  we  have 

THEOREi  3.1:  If  F  is  IFR  (DFR)  with  q^1  quantile  and  G 

th 

is  exponential  with  a  quantile  £  ,  then 

VSr„iC'(9lVi(i)  PG(§rn  (3'2) 


£ 


q 


(2r) 

£n(l-q) 


1 


when 
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Inequality  (3.2)  follows  from  Theorem  4.2  of  BP  (1966a).  (See  also 

Theorem  2.3  of  BP  (1966b).)  Note  that  the  coefficient  of  is 

greater  than  one  when  q  is  sufficiently  small;  i.e., 

2 

q  1  -  exp[-xa(2r)/2n] ,  so  that  the  consumer  is  protected  against 
accepting  material  not  meeting  the  quantile  life  goal. 

If  q  is  sufficiently  large  and  decision  rule  (3.1)  is  used, 
then  the  producer  rather  than  the  consumer  will  be  protected  when  the 
distribution  is  actually  IFRA.  More  precisely, 

THEOREM  3.2:  If  F  is  IFRA  (DFRA)  with  q  quantile  £  and 

L 

G  is  exponential  with  q  quantile  E,  3  then 

P-_ { 0  >  c’u°)k  }  >  (<)  p„{e  >  c'(£°)k  }  (3.3) 

F  rn  —  q  '  q  —  —  G  rn  —  q  1  q 


when 

; 

q - 2(n-r+l)£,n(l-q) 

Inequality  (3.3)  follows  from  Theorem  3.4  of  BP  (1966a). 

(See  also  Theorem  3.3  of  BP  (1966b).)  The  coefficient  of  £°  is 

q 

2 

less  than  one  for  large  q;  i.e.,  q  >_  1  -  exp[~xa(2r)/2(n-r+l)  ] . 
Hence  for  large  values  of  q,  the  producer’s  risk  is  controlled  but 


not  the  consumer's  risk. 


11 


The  analysis  above  points  up  the  fact  that  one  must  be  very 
careful  in  using  exponential  life  tests  based  on  censored  samples. 
Depending  on  the  objective  chosen,  either  the  producer  or  the 
consumer,  but  certainly  not  both  simultaneously,  can  be  protected 
if  the  distribution  is  actually  IFR  or  DFR. 


4.  SAMPLING  WITH  REPLACEMENT 

Another  device  to  save  time  spent  on  experimentation  is  to 
replace  failed  items  at  the  time  of  failure  with  new  items  until  the 
total  number  of  failed  items  reaches  a  fixed  number,  say  r.  Except 
for  this  restriction,  the  sampling  plan  is  the  same  as  in  Sections  2 
and  3.  We  shall  see  that  the  remarks  concerning  producer's  and 
consumer's  risks  are  even  more  applicable  in  this  case. 

Let  denote  the  time  of  the  i  failure  when  failed  items 

are  replaced.  The  maximum  likelihood  estimate  for  9  based  on  the 
exponential  assumption  is,  in  this  case 


9 

rn 


-  [nX*  +  n(X?  +  n(X?  -  X*)  +•••+  n(X*  -  X*  .)] 
r  1  1  2  1  r  r-1 


(4.1) 


the  total  time  on  test  divided  by  the  number  of  failures  observed 
(Epstein,  (1960b). 
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It  will  be  convenient  to  introduce  the  following  fictitious 
replacement  policy: 

Policy  A:  Replace  a  failed  item  with  a  good  item  of  the  same 
"age". 

th 

Let  denote  the  time  of  the  i  failure  under  this  policy, 

and 


*** 

0 

rn 


** 

"  Vl> 


/r, 


(4.2) 


the  total  time  on  test  divided  by  the  number  of  failures  observed. 
LEMMA  4.1:  If  F  is  IFR  (DFR) ,  then 


P  {0  >  x}  <  (>)  P  {§  >  x}. 

r  rn  —  —  —  r  rn  — 


(4.3) 


Proof.  Assume  F  is  IFR.  The  proof  is  by  induction  on  r. 

Note  that  (4.3)  is  obviously  true  for  r  =  1  and  all  n  >_  1.  Assume 

(4.3)  is  true  for  r-1  and  all  n  >  r-1.  Let  F,  denote  the 

—  In 

distribution  of  X,  , 

In 

Fu(x)  =  [F(x+u)  -  F(u)]/F(u) 

(the  distribution  of  an  item  of  age  u) ,  and 


T 

rn 


nX 


In 


n(X2n-Xln> 


+•••+  n(X  -X  ). 

rn  r-i,n 
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I 


Then 


w 

PjT  >  x}  =  /  V _{nX.  +•••+  (n-r+1)  (X  -X  .  )  >  x|nX,  =  uldF.  (7) 

F  rn  —  F  In  rn  r-l,n  -  '  In  In  n 


/ 


0  I>Fu«n-1>Xl,„-l  +-"+  ((-n-tr-l)  +  i)(xr.lin.i-xr.2,n.i)i»Wta(;) 


n 


/  P_  { T  ..  >  x-u}dF.  (-) 

F  r-l,n-l  —  In  V 


0  u 
n 


-  L  "f 


0  u 
n 


The  last  inequality  follows  from  Theorem  2.2. 


Also 


00 

itit  f  **  **  **  **  1  *  ,  U 

P{T  >  x}  =  /  P„{n(X_  -X.  )  +•••+  n(X  -X  .  )  >x-u  nX.  =  u}dF.  (-) 

rn  —  0  F  2n  ln  rn  r“^*n  ~  In  In  n 

OO 

«  /  PF  {T**  >  x-u}dF (~). 

r  r-l,n  —  in  n 

0  _u 

n 

By  the  induction  assumption 


it  it 

P_  (T  .  >  x-u}  <  P_  {T  .  >  x-u}. 

F  r-l,n  —  —  F  r-l.n  — 

u  u 


n 


n 


Hence 


00  00 

J  P_  {T  .  .  >  x-u}dF.  (-)  <  /  P_  {T**  >  x-u}dF-  (- 

n  F  r-1  ,n-l  -  In  n  -  Jn  F  r-l,n  -  lnvn' 


0  « 
n 


0  u 
n 


which  proves  (4.3). 


An  analogous  proof  holds  In  the  DFR  case. 


a  ** 

We  may  now  show  that  the  usual  estimate  0  under  censored 
J  rn 

sampling  with  replacement  is  stochastically  larger  than  the  corresponding 
estimate  0  under  censored  sampling  without  replacement. 

THEOREM  4.2:  If  F  is  IFR  (DFR),  then 


P  {0  >  x}  <  (>)  Pp{0  >  x}. 

r  rn  —  —  —  r  rn  — 


Proof .  Assume  F  is  IFR.  By  Lemma  4.1  it  is  sufficient  to 
show  that 


«**  a* 

P_{0  >  x}  <  P_{0  >  x}. 

F  rn  —  —  F  rn  — 


(4.4) 


To  show  this,  consider  the  set  of  n  tenewal  processes  generated  by 

it  AA 

successive  failures  and  their  replacements.  Let  N^(t)  (N^  (t)) 

denote  the  number  of  replacements  in  the  i^  renewal  process  under  the 
original  replacement  policy  (Policy  A).  Thus  to  show  (4.4)  we  need 
only  show 


&&  £ 
PK(Trnix)  <  PF(Trn>  x) 


or,  equivalently 
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0 


Since  life  lengths  in  different  renewal  processes  are  independent, 
we  need  only  show 


'fc'k  jl? 

VNi  (*)  1  r}  <  PpfN^x)  <  r}. 


(4.5) 


Let  { U  }  denote  a  renewal  process  with  underlying  distribution 
1  i=l 

00 

F,  IFR.  Let  {V  }  denote  a  sequence  of  random  variables  where 

1  i=l 

has  distribution  F,  and  given  ^  has  distribution 


FV  |V  ^ 
i1  i-1 


F(t+Vi_i>  -  F^) 


Then  (4.5)  is  equivalent  to 


P_,{N.  (x)  >  r}  =  P_,{V..  +•••+  V  <  x}  >  P_{U  +•  •  •+  U  <  x}  =  P  {N  (x)  >  r} 

r  1  —  r  1  r  —  —  r  1  r  —  r  1  — 

(4.6) 

To  prove  (4.6)  we  use  indue  .1  on  r.  Clearly  (4.6)  is  true  for 
r  =  1.  Suppose  it  is  true  tor  r-1.  Then 

X 

P{V  +•••+  V  <  x}  =  /  P_  {V  +• •  •+  V  .  <  x-u}dF(u) 
i  r  —  r  l  r-i 

0  u 

x  x 

>  /  P  {U  +•  •  •+  U  ,  <  x-u)dF(u)  >  /  P^UL  +•••+  U  ,  <  x-u}dF(u). 

— X  r  i  r-i  —  —  -v.fi  r-i  — 

0  u  0 

The  first  inequality  follows  from  the  induction  assumption  and  the 
second  from  the  IFR  property  that 
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F  (t)  >  F(t) 
u  — 

for  all  u  >  0  and  all  t  0. 

An  analogous  proof  holds  in  the  DFR  case.  |  | 

yif  ^ 

An  obvious  consequence  of  Theorem  4.2  is  that  E6rn  (f)  ® 

in  the  IFR  (DFR)  case. 


5.  TRUNCATED  LIFE  TEST  PLANS 

A  common  practice  in  life  testing  is  to  truncate  the  experiment 
at  a  pre-assigned  time,  say  T,  and  note  the  number  of  failures.  If 
the  number  of  failures  is  less  than  or  equal  to  an  acceptance  number  c, 
the  life  length  goal  is  considered  to  have  been  demonstrated;  otherwise, 
not.  Such  sampling  plans  have  been  discussed  by  Sobel  and  Tischendorf 
(1959)  for  the  exponential  case.  In  Barlow  and  Gupta  (1966)  such 
sampling  plans  are  considered  for  IFR  (DFR)  distributions  using  known 
bounds  for  these  distributions.  It  is  pointed  out  that  in  the  IFR  case 
we  have  to  test  beyond  the  mean  (or  quantile)  life  goal  to  protect  the 
consumer,  whereas  we  can  protect  the  producer  by  testing  for  a  time  less 
than  the  mean  (or  quantile)  life  goal.  Sampling  experiments  by  Zelen 
and  Dannemiller  (1961)  show  that  the  exponential  procedures  will  perform 
poorly  when  the  underlying  distribution  is  Weibull  if  the  truncation 
time,  T,  is  less  than  the  mean  life  goal. 


Bartholomew  (1963)  considers  truncated  life  test  plans  for 
the  exponential  distribution  where  the  times  of  failures  occurring 
before  T,  the  truncation  time,  are  assumed  known.  Assuming  at 
least  one  failure  has  occurred  by  time  T,  the  maximum  likelihood 
estimate  for  the  mean,  6,  is 

eCO  *  7  E  fT  -  a.tt-Y  )] 
r  i  11 

where  Y,  <  Y_  <  •••  <  Y  denotes  an  ordered  sample  from  an 
1  —  Z  —  —  n 

exponential  distribution  and 

if  Y,  <  T 
i  — 

otherwise, 

n 

and  r  =  E  a.j  •  The  distribution  of  this  statistic  is  computed  in 
1 

Bartholomew  (1963).  As  in  the  censorship  case,  an  acceptance  test 
based  on  this  statistic  will  favor  the  producer  if  the  true  distribution 
is  IFR  and  T  is  not  too  large. 

THEOREM  5.1:  If  F  is  IFR  (DFR)  with  mean  6,  T  <  6,  and 
G(t)  =  1  -  e  then 

Pp{ 6 (T)  >_  x | r  >  1}  1  (±)  PG{§(T)  >  x |r  >  1} 


for  all  x  >  0. 
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Proof .  Assume  F  is  IFR.  Let  *  -0  log  F(X^)  and  let 

1  if  Yi<T 

bi  * 

0  otherwise. 

Since  G(x)  crosses  F(x)  to  the  right  of  0  and  from  above 

(Barlow  and  Proschan  (1965),  Chapter  2),  we  see  that  <  T  implies 

Y  <  T,  and  hence  a±  <  Thus  if  a1  =  1,  then  bi  *  1  and  so 

n 

T-X.  <_  T-Y . .  Hence  if  £  a.  1, 
ll  i  i 

n  n 

£  [T-a1(T-X1>]  £  [T-bi(T-Yi> ] 

n  —  n 

Ea,  Eb, 

1  1  1 

By  Lemma  1,  page  73,  Lehman  (1959),  it  follows  that 

Pp(0(T)  _>  x|r  >_  1}  >_  Pg{0(T)  >_  x|r  >_  1} 

for  all  x  ^  0  and  T  <  0. 

A  similar  proof  holds  in  the  DFR  case.|| 

Bartholomew  (1957)  uses  the  following  estimate  for  0  in  the 


exponential  case: 
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e(T)  = 


,  r 

;Ey.  +  (n-r)T 
1  1 

nT 


if  r  >  0 

if  r  =  0. 


Note  that  for  r  >  0,  this  is  the  maximum  likelihood  estimate 
discussed  above.  Then 


E[e(T)]  -  e  -  cov^lq»  >  e 
1-exp  (-  -g-) 

since  r  and  0(T)  are  negatively  correlated,  Bartholomew  (1957). 

Thus  the  estimate  is  positively  biased  when  the  underlying  distribution 
is  exponential.  In  BP  (1966a)  it  is  shown  that 

r  jr 

Ep(i;xi  +  (n-r)T]  >  (<)  Eg  [  ^  Y±  +  (n-r)T] 

when  F  is  IFRA  (DFRA) ;  T  is  not  required  to  be  <0  as  in 
Theorem  5.1.  This  would  suggest  that  the  estimate  is  even  further 
biased  when  the  underlying  distribution  is  IFRA. 

Inverse  Binomial  Sampling.  Nadler  (1960)  has  considered  the 
following  type  of  sampling:  An  item  having  life  distribution  F  with 
mean  0  is  put  on  test  until  it  fails  or  time  T  has  elapsed;  at  this 
time  the  item  is  replaced  by  a  fresh  item.  This  is  repeated  sequentially 
until  r  actual  failures  are  observed.  The  number,  N,  of  items  that 
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have  to  be  tested  until  r  actual  failures  are  obtained  is  a  random 

variable.  Let  X, ,X„,...,X  denote  a  sequence  of  independent 

12  n 

random  variables  with  distribution  F.  Let 


'  X  if  <  T 

Z.  =  • 

1 

,  T  otherwise, 


so  that 
when  F(x)  * 


is  the  time  on  test  of  the  i*"*1  item.  Nadler  showed  that 

-X  /  0 

1-e  ,  a  sufficient,  unbiased  estimate  of  0  is 


er(T) 


l 

r 


£ 


1 


L*  ,  • 
i 


A 

We  compare  6  (T)  with  related  statistics  in 

THEOREM  5.1:  If  F  ir  IFR  (DFR)  with  mean  0,  then 


pF{er(T)  >_  x}  >  (<)  PF{err  >  x}  (5.1) 

for  all  x  >_  0,  and  hence 

EF[§r(T)]  >  (<)  6. 

Ifj  in  addition }  T  <  6,  then 

Pr  (•'('!)  >  x}  >  (<)  P  { G  (T)  >  x}  =  P_{0  >  x) 

F  i  —  —  —  G  r  —  G  rr  — 


(5.2) 
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for  all  x  >_  0. 

Proof .  Assume  F  is  IFR.  Let  denote  the  time  between 

the  (i-l)st  and  i^  actual  failure  following  an  inverse  binomial 

1  r 

sampling  policy.  Note  that  -  =  6  (T) .  Then  for 

r  1  r 

nT  <_  x  <  (n+l)T, 

P[Wt  >x]  =  [ F (T)  ]n[F(x-nT)  ]  >_  F(x)  =  P(X±  ^  x] ; 
the  inequality  follows  from  the  fact  that  F  is  IFR.  Hence 


proving  (*>.1). 

To  show  (5.2),  let  Yj  =  -log  F(Xi>.  Note  that  Yi  <_  if 

X.  <  T  <  0.  It  follows  that  if  N*  denotes  the  number  of  Y.  's 
x  —  i 

must  be  examined  to  obtain  r  values  each  less  than  T,  we  see 
that  N'  <_  N.  Therefore,  for  each  sample, 

N'  N 

£  z;  <  Ezi. 

i  i 

where 

Yi  if  Y  <_  T 

T  otherwise. 


which 
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(5.2)  follows  from  Lemma  1,  page  73  of  Lehman  (1959).  | 


6.  BOUNDS  ON  ESTIMATES  FOR  THE  RELIABILITY  FUNCTION 

In  this  section  we  obtain  a  lower  bound  for  the  expected  value  of 
two  estimates  of  reliability  appropriate  under  the  exponential  assumption, 
when  the  distribution  is  actually  IFR.  Assume  that  a  censored  sampling 
plan  is  followed,  as  in  Section  2. 

The  minimum  variance  unbiased  estimate  for  R(t)  *  F(t)  (t  is  fixed) 
under  the  exponential  assumption  is 


Rx(t) 


0 


if  t  <  Z 

r 

otherwise, 


r 

where  Z  =  ^  X  +  (n-r)X  ,  the  total  time  on  test.  For  a  discussion 
r  1  1  r 

of  such  minimum  variance  unbiased  estimates,  see  Tate  (1959).  Then,  under 
the  exponential  assumption  with  the  mean  0  taken  to  be  1  for 
convenience,  Z  has  density 


r-1  -y 

«.w  ’  fcifr 


We  obtain  a  lower  bound  for  R^(t)  in 
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THEOREM  6.1:  If  F  is  IFR  with  ean  6-1  and  t  <  0,  then 


E[R  (t)  >_  f  I  1  -  £  I  g  (y)dy  +  (l-t)r  1  /  g  (y)dy. 

1  t  1  y  1  1 


Proof .  Let  <J>  ^(x)  =  -log  F(x).  Then  <J>  is  concave,  increasing 
and  <j)(0)  =  0.  If  is  the  i^  order  statistic  in  a  sample  of  n 

from  an  exponential  distribution  with  mean  1,  then  =  <t(Y^)  is 
distributed  as  the  ith  order  statistic  from  F.  Furthermore, 


Z  =  E<j.(Y  )  +  (n-r)  <J> (Y  ) 
1 


r  r  1  * 

1  *  I  EYi  +  (n-r)Yr  j  =  d»(Zr). 


Therefore , 


Rx(t)  -  1 


.  r-1 


1  -  TT~ 


if  t  <  Z 


otherwise 


1  Rx(t) 


-  < 


1  -  t/4.(Zr) 


r-l 


if  t  <  <j>(Zr) 


otherwise . 


Since 
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<P(y) 


y 

1 


y  <  i 
y  1 1. 


we  have 


E[R^(t) ]  > 


r-1 

gr(y)dy  + 


/  (l-t)r-1g  (y)dy.  |  | 
1 


The  maximum  likelihood  estimator  for  R(t)  under  the  exponential 
assumption  is 


R2(t)  =  exp[-t/Grn] 

A 

where  0  is  defined  in  Section  2.  Pugh  (1963)  has  shown  that  under  the 

A 

exponential  assumption  R2(t)  is  negatively  biased  when  the  true 
reliability  R(t)  >  ^  ~  .368.  Assuming  F  is  IFR  we  can  obtain  a  lower 
bound  on  (t)  ] . 

THEOREM  6.2:  If  F  is  IFR  with  mean  0  =  1,  then 

1  tr  oo 

-  /  e  y  +  e  tr  /  gr(y)dy. 

0  1  r 


The  proof  parallels  that  of  Theorem  6.1. 
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